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1.0  INTRODUCTION 


Recently,  a  study  group  sponsored  by  the  Amiv  Research  Office  developed  and  published 
an  overall  basic  research  plan  to  guide  advanced  research  on  key  nitramine  propellant  ignition  and 
combustion  issues  which  must  be  addressed  to  allow  the  systematic  exploitation  of  this  exciting 
class  of  energetic  materials. 1  Subsequent  workshops  organized  by  the  Office  of  Naval  Research 
and  the  JANNAF  Combustion  Subcommittee  on  Nitramine  Propellants  have  further  refined  the 
outstanding  issues.  The  work  reported  here  is  a  response  to  a  need  identified  by  all  these  groups, 
the  development  of  new  diagnostic  techniques  to  monitor  the  progress  of  decomposition  and 
ignition  kinetics  in  the  condensed  phase. 

There  have  been  very  few  direct  measurements  of  condensed  phase  decomposition  in 
burning  gun  propellants.  Our  technique,  motivated  by  recent  advances  in  electro-optical  material, 
is  the  use  of  embedded  optical  fibers.  We  have  begun  by  investigating  the  use  of 
infrared-transmitting  fibers  in  a  direct  absorption  experiment,  in  which  a  gap  between  a  pair  of 
embedded  fibers  may  allow  spectroscopic  detection  of  decomposition  products  within  a  strand  of 
solid  propellant. 

Although  ours  is  the  first  spectroscopic  investigation  in  the  solid  phase  of  burning 
propellants,  it  has  benefitted  from  a  number  of  related  research  programs  (a  recent  comprehensive 
review  of  experimental  and  theoretical  work  is  found  in  Reference  2).  Infrared  absorption  has 
been  used  to  study  the  decomposition  of  propellant  ingredients,  including  the  nitramines  HMX  and 
RDX  V7  as  well  as  other  materials.8  These  studies  all  used  Fourier  transform  infrared 
spectrometers  to  measure  the  changing  absorption  of  thin  films  of  the  materials  as  they  underwent 
rapid  heating^'6  or  photolysis  by  an  ultraviolet  lamp.7  In  particular,  the  spectra  reported  in 
Reference  3.  in  which  the  absorption  bands  of  NiO  and  CO2  can  be  seen  as  these  gases  accumulate 
in  bubbles  forming  in  confined  thin  films  of  HMX  and  RDX.  did  much  to  encourage  us  in  our 
approach.  Other  programs  providing  useful  information  have  included  the  work  by  Miller  on 
embedding  line-wire  thermocouples  in  burning  propellant  strands,9  the  ongoing  work  at  Sandia 
National  Laboratories  involving  infrared  spectra  o I  heated  RDX  films,10  and  work  by  Lit/.inger 
and  co-workers  on  mass  spectrometer  sampling  of  laser-heated  propellant  samples.1 1 


2D  EXPERIMENTAL  TECHNIQUES 


Our  goal  was  to  develop  a  technique  of  last-response  spectral  infrared  fiber  probing  of 
burning  propellant  strands,  and  demonstrate  it  in  mechanism  studies  of  condensed  phase  processes 
m  nitraminc  propellant  combustion.  In  early  stages  of  this  work,  we  focused  on  observing 
chemical  changes  in  solid  propellant  by  infrared  absorption,  embedding  a  pair  of  fibers  separated 
by  a  sliver  of  propellant.  We  will  not  report  on  that  work  here,  but  instead  confine  our  discussion 
to  a  second  experimental  configuration  which  has  been  much  more  successful.  In  these 
experiments,  the  gap  between  fibers  (inside  the  propellant  strand)  was  lengthened  to  0.2  to  0.5  cm 
and  left  empty.  This  air-filled  cavity  or  observation  volume  could  then  fill  with  gas-phase 
decomposition  products  as  the  flame  front  approached.  Observations  including  when,  how  fast, 
and  to  what  levels  it  was  filled  could  be  used  to  support  hypotheses  as  to  the  source  of  the 
decomposition  gases. 

Previous  studies. I‘^’ 1  indicated  that  nitraminc  propellant  ingredients  such  as  RDX  may 
decompose  into  a  number  of  gaseous  species  including  NyO.  H2CO.  NOy.  HCN.  HyO.  NO.  and 
COy.  The  fiber  length  used  here  of  about  1  m  means  the  transmission  band  of  the  fiber  matches 
the  response  curve  of  the  detectors,  both  extending  from  less  than  2  to  about  5  pm.  This  allows 
sensitive  detection  of  NyO  in  the  Vy  band  around  4.7  pm.  as  well  as  COy  in  the  band  around 
4.3  pm,  HyCO  around  3.5  pm  and.  with  a  laser  source,  NO  around  5.2  pm.  Of  these.  NyO  was 
chosen  as  the  first  example  because,  of  those  species  expected  to  give  strong  infrared  absorption 
signals  with  the  first  version  of  our  apparatus,  it  was  of  greatest  interest  for  the  definition  of 
chemical  mechanisms  of  decomposition.  Extensions  of  optical  fiber  probe  techniques  to  other 
visible  and  infrared  wavelengths  to  quantify  other  potential  decomposition  species  including  NOy 
are  being  explored  in  a  follow-on  effort. 

A  schematic  diagram  of  the  experiment  reported  here  is  illustrated  in  Figure  1.  The 
chopping  wheel  modulates  the  light  from  an  infrared  source  (a  globar).  and  so  the  transmitted 
signal,  at  approximately  2(H)()  H/„  The  light  is  focused  onto  a  fluoride  glass  optical  fiber  (Infrared 
Fiber  Systems)  of  200  pm  nominal  diameter.  The  propellant  strand  is  mounted  inside  a  high 
pressure  strand  burner  (shown  in  Figure  2).  The  embedded  optical  fibers  pass  through  the 
baseplate  of  the  strand  burner,  as  do  leads  connecting  a  variac  to  a  V-shaped  segment  of  nichrome 
wire  pressed  against  the  top  of  the  strand  to  provide  ignition.  In  the  observations  reported  here, 
the  gaps  between  the  source  and  detector  fibers  are  in  the  0.2  to  0.5  cm  range.  Light  collected  by 
the  liber  on  the  detector  side  is  divided  by  a  50/50  beamsplitter  and  passes  through  a  pair  of  narrow 
band  infrared  interference  filters  before  being  focused  onto  InSb  infrared  detectors. 
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Figure  I.  Schematic  Drawing  of  Infrared  Fiber  Optic  Absorption  Experiment. 

A  data  acquisition  computer  reads  a  signal  point  from  each  infrared  detector  at  each  on  and 
off  point  of  the  chopper.  Computer  data  acquisition  for  each  chopper  cycle  is  triggered  by  the 
synchronization  signal  from  the  commercial  chopper  controller.  The  signals  from  the  infrared 
detectors  and  pre-amplifiers  pass  through  buffer  amplifiers  to  minimize  the  effects  of  the  analog  to 
digital  conversion.  This  also  allows  optimization  of  the  time  response  of  the  signal,  so  the  single 
data  point  the  computer  reads  is  representative  of  the  entire  half-cycle  of  the  chopper. 

The  strand  burner  shown  in  Figure  2  is  a  copy  of  a  device  in  use  at  the  U.S.  Army 
Research  Laboratory,  Aberdeen  Proving  Ground,  with  some  modifications  to  accommodate  fiber 
optics.  The  only  propellant  studied  to  date  has  been  XM39.  composed  of  15(7<  RDX  and  25 (7> 
binder.  Its  detailed  composition  is  shown  in  Table  1 .  Experiments  were  performed  at  2  pressures. 
I  and  6  atm.  Although  both  pressures  are  far  from  those  encountered  in  rocket  combustion 
chambers.  let  alone  in  cannon  shells,  it  turns  out  that  the  variation  in  burning  rate  between  the  two 
pressure  regimes  is  large  enough  to  give  qualitatively  different  observations  and  hence  useful 
information. 
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Fiuure  2.  Cross-Section  of  Strand  Burner. 


Table  1.  Composition  of  XM39  Propellant 


Component  Name 

Chemical  Formula 

Mass  Fraction 

Meltine  Point 

RDX  (5  micron) 

c,h6n6o6 

0.76 

190°C 

Cellulose  Acetate 
Butyrate  (CAB) 

C,sH2208 

0.12 

2 1 5-225  °C 

Nitrocellulose  (NC) 

ChH75505(N02)245 

0.04 

160-170  °C 

Acetyl  Triethyl 

Citrate  (ATC) 

(  |4E12208 

0.076 

(B.  P.  132  °C) 

Ethyl  Centralite  (EC) 

c17h20on2 

0.004 

72  °C 

XM39  will  not  bum  in  an  inert  atmosphere  at  1  atm.  so  we  used  secondary  combustion  of 
the  gaseous  decomposition  products  with  air  to  support  the  primary  combustion  of  the  solid 
propellant.  A  combination  of  nitrogen  addition  to  the  air  and  a  large  enough  How  velocity  around 
the  strand  was  used  to  prevent  burning  down  the  side  of  the  strand.  At  atmospheric  pressure,  we 
used  a  gas  flow  through  the  central  inlet  composed  of  235  STP  cm  Vs  air  and  40  STP  cm  Vs  N2. 
for  a  gas  velocity  around  the  strand  of  roughly  3.5  cm/s.  At  6  atm.  the  air  and  N2  Hows  were 
increased  to  2470  and  825  STP  cm  Vs,  respectively,  changing  the  velocity  of  the  gas  flow  around 
the  strand  to  around  7  cm/s.  Under  these  gas  How  conditions,  the  strand  burns  w'ilh  a  roughly  Hal 
melt  layer  and  a  slightly  concave  burning  surface.  From  measurements  discussed  below,  we 
found  propellant  burning  rates  of  about  0.015  cm/s  at  I  atm.  and  0.041  cm/s  at  6  atm. 

A  high-resolution  video  system  is  used  to  record  each  burn.  The  beginning  of  data 
collection  is  accompanied  by  a  computer  signal  to  the  video  recorder  which  starts  its  time  clock. 
With  a  second  lens,  the  video  camera  is  also  used  as  a  microscope  during  fiber  insertion. 

The  libers  are  inserted  between  two  sections  of  a  propellant  strands,  cut  with  a  diamond  saw. 
The  liber  insertion  procedure  begins  with  the  formation  of  a  groove  in  both  propellant  pieces.  The 
strand  segment  is  placed  in  a  hole  in  a  warmed  metal  block.  At  the  bottom  of  the  hole  is  a  wire 
stretched  across  its  diameter.  The  wire  thickness  is  chosen  so  that  heavy  pressure  on  the  strand 
leaves  the  proper  indentation  on  its  face.  Using  a  knife,  grooves  for  thermocouples  (0.0127  cm 
diameter,  chromel-alumel)  are  cut  to  within  approximately  0. 1  cm  of  the  fiber  groove.  Then,  after 
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I  ho  l  i  hers  arc  positioned  in  the  groove,  they  are  pressed  into  the  groove  using  the  top  propellant 
piece,  the  thermocouples  are  inserted,  and  both  libers  and  thermocouples  are  secured  inside  the 
strand  with  cyanoacrylate  adhesive  at  the  entry  points.  Figure  3  illustrates  the  technique.  The 
reassembled  strand  has  a  cylindrical  cavity  or  observation  volume  with  a  diameter  of  about  0.02  cm 
and  a  length  of  0.2  to  0.5  cm.  with  thermocouple  junctions  on  one  or  both  sides  ol  it.  a  lew  liber 
diameters  outside  the  wall  of  the  cavity  and  from  0.005  cm  to  0.02  cm  below  the  plane  containing 
the  liber. 


Figure  3.  Fiber  Insertion  Procedure. 


3.0  EXPERIMENTAL  ANALYSIS 


3. 1  Nitrous  Cxide  Quantification 

The  509?  transmission  points  of  the  filter  used  to  detect  NoO  were  at  about  2 1 85  and 
2238  cnr1.  while  those  of  the  reference  filter  were  2093  and  2196  cnr1.  Since  the  half-height 
points  of  the  NfO  absorption  band  are  at  around  2184  and  2254  cnr1,  it  is  clear  not  only  that  the 
NxO  filter  closely  overlaps  the  N2O  band,  but  also  that  the  reference  filter  detector  will  be 
somewhat  sensitive  to  the  presence  of  N2O  in  the  absorption  path.  Initially,  we  simply  included 
this  sensitivity  in  the  analysis  relating  the  ratio  of  the  two  filter  intensities  to  N2O  mole  fractions. 

In  later  experiments  a  short  cell  filled  with  one  atmosphere  of  N2O  was  placed  in  front  of  the 
reference  filter,  so  that  light  reaching  the  reference  detector  only  had  wavelengths  which  are  not 
absorbed  by  N2O,  thus  eliminating  the  sensitivity. 

The  filter  responses  were  calibrated  computationally,  using  Fourier  transform  infrared 
spectrometer  measurements  of  the  bandpass  filter  transmission  curves,  and  N2O  spectral 
absorbance  values  derived  from  the  HITRAN  absorption  line  database12.  Although  it  was  checked 
against  observations  made  by  filling  the  strand  burner  with  a  range  of  pressures  of  N2O  while  the 
gap  between  fibers  was  still  uncovered,  the  computational  calibration  is  the  primary  basis  for  the 
observed  mole  fractions  reported  here,  because  it  most  easily  allows  us  to  take  into  account 
changes  in  the  in-band  absorption  due  to  changes  in  line  strengths  and  widths  with  temperature. 

The  transmitted  intensity  signals  are  obtained  from  computer  files  of  the  detector  outputs  by 
differencing  each  pair  of  chopper-open  and  chopper-closed  values.  These  intensity  records  are 
then  averaged —  an  average  of  63  pairs  of  points  yields  a  30  Hz  data  record,  and  for  noisier  data 
10  Hz  or  5  Hz  averages  are  used.  After  the  fiber  has  burned  through,  the  remaining  data  stream  is 
averaged  to  give  a  zero  level  which  is  subtracted  from  each  detector  signal  before  the  ratio  is  taken. 
The  value  of  this  ratio,  relative  to  its  value  before  ignition,  gives  a  fractional  transmission  in  the 
N2O  band.  It  is  compared  with  the  calibrations  described  above  to  directly  yield  the  product  of  the 
pathlength  and  the  N2O  number  density  and.  with  further  assumptions,  the  N2O  mole  fraction. 

3.2  Flame  Front  Position  Measurement 

The  video  record  ol  each  burn  was  examined  to  yield  a  record  of  the  position  of  the  burning 
surface,  synchronized  with  the  computer  records  of  infrared  transmission  and  thermocouple 
temperatures.  Flame  front  positions  were  typically  measured  at  intervals  of  one  second  or  more 
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with  a  lime  uncertainty  of  less  than  the  frame  period  of  1/30  s  and  a  position  uncertainty  of  about 
0.02-0.03  cm  (the  raster  spacing  lor  the  range  of  image  magnifications  used).  These  position 
histories  typically  exhibit  quite  steady  burn  rates.  Some  burns  do  have  irregularities  near  the  fiber 
insertion  point.  These  irregularities  may  in  some  cases  be  due  to  a  change  in  the  bum  rate  across 
the  entire  face  of  the  propellant  strand,  but  many  apparent  irregularities  are  due  to  difficulties  in 
seeing  the  center  of  the  strand,  such  that  the  point  chosen  to  represent  the  burning  surface  is  on  a 
front  edge.  An  edge  region  may  bum  slower  or  faster  as  the  fiber  insertion  point  is  approached, 
because  of  the  heat  transfer  and  air  How  perturbations  of  the  fibers,  thermocouples,  adhesive,  and 
the  cut  through  the  strand. 

There  could  be.  therefore,  a  difference  of  as  much  as  0. 1  cm  in  some  cases  between  the 
position  in  the  transcription  of  the  video  record  and  the  actual  distance  from  the  flame  front  in  the 
center  of  the  strand  to  the  fiber  insertion  plane.  To  give  a  measure  of  how  much  these  differences 
affect  ensemble  averages  of  a  quantity  such  as  the  liber-to-llame-front  distance  w  hen  N2O  first 
appears  in  the  observation  volume,  we  will  report  this  distance  based  both  on  the  point-by-point 
video  record,  and  on  a  least  squares  fit  to  the  entire  position  history.  The  two  positions  agree  far 
from  the  fiber  insertion  point,  and  if  the  burning  rate  of  the  center  of  the  strand  were  more  regular 
than  that  of  the  edges,  we  might  expect  the  linear  fit  to  give  somewhat  more  reliable  fiber-to-flame- 
lront  distances. 

3.3  Propellant  Temperatures  Measurements 

As  noted  above,  the  thermocouples  had  0.0127  cm  wire  diameters,  and  the  diameter  of  the 
bead  at  the  junction  is  more  than  double  this  value.  However,  we  are  convinced  that  the  major 
source  of  unccilainty  in  the  thermocouple  readings  is  not  this  relatively  large  size  or  issues  of  good 
thermal  contact  with  the  propellant  or  heat  transfer  out  through  the  thermocouple .  but  rather  the 
t occasional)  irregularity  in  the  burning  surface  in  the  vicinity  of  the  fibers.  As  discussed  above, 
the  burning  surface  far  from  the  fibers  is  relatively  perpendicular  to  the  axis  of  the  strand  and  is 
slightly  concave.  As  the  flame  bums  through  the  cut  in  the  strand  and  around  the  adhesive  beads 
on  the  strand  surface  which  secure  the  fibers  and  thermocouples,  the  slope  of  the  burning  surface 
can  become  larger,  such  that  the  thermocouple  bead  could  be  in  a  warmer  or  colder  region  of 
propellant  than  the  fiber  gap  of  order  0. 1  cm  away  from  it.  It  should  also  be  said,  though,  that  in 
our  burns  in  which  two  thermocouples  were  inserted,  the  two  temperature  traces  agreed  well. 
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4.0  EXPERIMENTAL  RESULTS 


4. 1  Example  Data  Traces 

Figures  4  and  5  present  examples  of  the  data  traces  near  the  time  of  NSO  appearance  for  the 
two  pressure  regimes  studied  ( 1  and  6  atm).  The  figures  show  the  two  detector  intensities  (both 
normalized  io  one  at  their  maximum  values),  their  ratio,  the  embedded  thermocouple  temperature, 
and  the  two  measures  of  distance  between  the  burning  surface  and  the  fiber.  In  both  cases,  the 
N2O  (as  indicated  by  the  drop  in  the  ratio  of  detector  intensities)  appears  and  reaches  a  fairly 
constant  level  in  a  time  between  0. 1  and  0.2  s.  It  must  be  pointed  out  that  the  two  pressures  differ 
in  how  often  this  abrupt  N2O  appearance  is  observed,  since  it  is  seen  in  5  of  6  cases  at  6  atm. 
while  in  5  of  7  1  atm  cases  the  N2O  appears  over  times  ranging  from  0.3  to  2.0  s. 

If  we  assume  that  the  NoO  is  at  the  temperature  of  the  embedded  thermocouples  (as 
mentioned  above,  an  assumption  with  some  uncertainty,  but  we  have  no  better  source  of  gas 
temperature  information)  we  can  use  the  known  absorption  path  lengths  to  calculate  NoO  mole 
fractions,  which  are  0.05  for  the  6  atm  case  in  Figure  4  and  0. 1 8  for  the  1  atm  case  in  Figure  5.  If 
we  make  the  conservative  assumption  that  the  ranges  of  thermocouple  temperatures  at  which  N20 
appears  (discussed  below)  represent  the  uncertainty  in  the  temperature  of  the  N20  in  the 
observation  volume,  then  from  our  precalculated  filler  calibration  factors  we  find  that  at  1 
atmosphere  this  translates  into  a  relative  uncertainty  of  less  than  ±0.3  in  N20  mole  fraction,  while 
at  6  atm.  the  relative  uncertainties  in  N2O  vary  within  the  range  of  ±  0.05  to  ±  0.20  as  the 
transmission  varies  from  0.8  to  0.2. 

After  its  initial  abrupt  appearance,  the  N2O  level  (and  the  infrared  transmission  through  the 
observation  volume)  is  maintained  for  an  average  of  about  0.4  s  in  the  6  atm  cases  and  for  between 
1  and  2  s  in  the  1  atm  cases  (in  the  1  atm  cases  with  more  gradual  N2O  appearance  the  time 
between  first  appearance  and  loss  of  transmission  is  also  on  the  order  of  1  s).  We  will  report  two 
time  segments  for  the  abrupt  appearance  cases,  referring  to  them  as  fill  times  and  steady  state 
times,  in  spite  of  the  facts  that  in  other  cases  they  are  pari  of  a  continuous  event  of  N2O  appearance 
in  the  observation  volume,  and  that  N2O  absorption  after  the  initial  fill  is  only  relatively  constant. 
On  this  latter  point,  we  can  note  that  in  Figure  4  a  change  in  N2O  temperature  from  1(X)  to  140  °C 
(following  the  rising  thermocouple  reading)  would  be  expected  to  cause  the  transmission  to 
increase  from  0.74  to  0.78 —  in  fact,  if  anything  it  decreases,  suggesting  that  the  NSO  mole 
fraction  is  increasing  in  the  observation  volume,  although  at  a  much  slower  rate. 


Alter  some  time  at  a  tairly  steady  NAO  level,  a  second  relatively  abrupt  event  brings  the 
transmission  to  zero.  It  can  be  seen  in  this  region  that  the  absolute  transmission  can  drop  by  an 
order  of  magnitude  while  the  ratio  of  the  two  transmissions  remains  constant.  Therefore,  die 
experiment  is  concluded  by  events  whose  effects  on  the  transmission  have  no  spectral  dependence 
over  the  range  of  the  bandpass  filters,  such  as  misalignment  of  the  fibers  or  entry  of  liquid  into  the 
observation  path. 

Figures  4  and  5  exemplify  a  second  difference  between  the  two  pressure  regimes,  that  the 
appearance  of  the  N2O  takes  place  in  different  temperature  ranges  —  at  under  1(K)  °C  in  the  6  atm 
example,  but  when  the  thermocouple  reading  is  about  250  °C  in  die  1  atm  case.  Finally,  the  values 
of  fiber  insertion  plane  to  burning  surface  distance  at  NyO  appearance  seen  in  the  two  figures  are  in 
fact  representative  of  the  range  seen  in  all  cases.  At  6  atm.  interpolation  in  the  point-by-point 
record  gives  a  distance  of  0.12  cm  while  using  the  linear  fit  through  all  position  points  gives  a 
distance  of  0. 1 3  cm.  At  1  atm.  the  two  distance  values  are  0.2 1  and  0.24  cm. 

As  mentioned  above.  Figure  4  is  representative  of  the  majority  of  the  6  atm  cases,  but 
Figure  5  is  not  representative  of  most  of  the  1  atm  cases,  which  instead  look  like  the  example  in 
Figure  6.  Here,  the  NyO  appears  gradually,  over  a  umc  period  of  at  least  a  second,  before 
transmission  is  lost.  The  temperature  and  flamc-to-fiber  distance  ranges  where  transmission  is  lost 
are  not  significantly  different  from  those  seen  in  the  abrupt  cases. 

The  1  atm  cases,  with  a  much  smaller  change  in  transmission  for  the  same  N2O  mole 
fraction,  usually  had  a  signal-to-noise  level  like  that  shown  in  Figure  6.  Without  fitting  a  smooth 
curve  to  the  observed  transmission  it  is  difficult  to  see  that  there  is  a  systematic  decrease  in 
transmission  corresponding  to  NAO  appearance.  The  one  1  atm  case  in  which  NAO  appearance  is 
easily  visible,  shown  in  Figure  7.  is  also  one  of  two  cases  to  show  a  second  phenomenon:  the 
subsequent  disappearance  of  N2O.  again  on  the  relatively  slow  time  scale  of  a  few  seconds. 

This  phenomenon  was  also  seen  in  the  6  atm  case  shown  in  Figure  X.  Once  again,  tne  N2O 
both  appears  and  disappears  on  a  relatively  slow  time  scale.  Furthermore,  a  second,  more  abrupt 
and  higher  concentration  appearance  of  N2O  may  occur  just  before  transmission  loss.  (Note  that 
the  data  set  of  Figure  X  is  clearly  one  of  those  taken  before  the  addition  of  the  NAO  gas  filler  in 
front  of  the  reference  detector,  since  both  detector  signals  respond  to  the  appearance  of  N2O.)  The 
cases  in  Figures  7  and  X  share  two  characteristics  which  may  help  to  explain  why  only  they  show 
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- Ratio  of  Detector  Intensities  (N20/Reference) 

- Normalized  N20  Detector  Intensity 

- Normalized  Reference  Detector  Intensity 

- Embedded  Thermocouple  Reading,  Degrees  C 

—  —  Point-by-Point  Flame  to  Fiber  Distance,  cm 

—  •  —  Linear  Fit  to  Flame  to  Fiber  Distance,  cm 


Figure  4.  Example  Data  for  XM39  Burning  at  6  atm  Pressure. 

NiO  disappearance:  they  both  maintained  good  infrared  transmission  for  several  seconds  after  NSO 
appearance,  and  their  thermocouple  temperature  readings  in  the  steady  state  region  were  the  highest 
in  each  pressure  set. 

In  other  experiments,  we  found  that  the  fluoride  glass  fibers  can  maintain  high 
transmission  even  when  heated  well  above  the  softening  point.  The  loss  of  transmission  depends 
on  mechanical  deformation,  not  on  the  fiber  temperature  per  se.  In  the  strand  experiments,  another 
factor  would  be  misalignment  of  the  two  fibers  once  the  surrounding  material  has  melted.  Thus, 
with  improved  support  for  the  fiber  it  might  be  possible  to  obtain  observations  such  as  those  in 
Figure  7  and  X  routinely. 

The  final  example,  in  Figure  9,  was  done  to  add  support  to  the  hypothesis  that  abruptly 
appearing  NoO  is  not  evolved  from  the  walls  of  the  observation  volume  but  arrives  from  regions 
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Figure  5.  Example  Data  for  XM39  Burning  at  1  atm  Pressure.  (Same  key  as  Figure  4).  In  this 
run,  an  N20  gas  filter  cell  was  not  used,  so  both  detectors  respond  to  N20. 
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Figure  6.  Example  Observations  for  XM39  Burning  at  1  atm  Pressure,  Showing  Slower  N->0 
Evolution.  (Same  key  as  Figure  4.) 
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19  Burning  at  1  atm  Pressure.  Showing  Slow  Appearance  and  Disappt 
.  (Same  key  as  Figure  4.) 
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Figure  X.  XM39  Burning  at  6  atm  Pressure,  Showing  Slow  Appearance  and  Disappearance  of 
N2O.  (Same  key  as  Figure  4.) 
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Figurc  9.  XM39  Burning  at  6  atm  Pressure,  With  Kapton  Film  Inserted  Above  Fiber  Gap. 

closer  to  the  flame  front  through  suddenly  opening  cracks.  In  two  cases,  experiments  were  set  up 
exactly  like  the  scries  of  6  atm  experiments  which  yielded  sudden  N20  appearance,  except  that 
before  the  two  segments  of  propellant  strand  were  joined,  a  0.0025  cm  thick  sheet  of  Kapton  (a 
strong,  high-temperature  plastic)  was  placed  on  top  of  the  fibers  (covering  about  the  central  half  of 
the  strand  surface).  With  this  more  solid  "roof",  no  abrupt  appearance  of  N20  was  observed. 
Instead,  the  N20  effused  slowly  around  the  plastic  sheet  and  into  the  observation  volume. 

4.2  Numerical  Results 

Table  2  lists  averages  or  ranges  of  some  of  the  parameters  which  can  be  derived  from  the 
data  traces  exemplified  by  Figures  4  through  8.  Beginning  with  the  bum  rates  derived  from  the 
video  records,  it  goes  on  to  address  five  lime  points  in  the  history  of  a  bum:  the  appearance  of  N20 
in  the  observation  region,  then  two  points  at  which  the  temperature  history  changes  its  functional 
form,  the  point  at  which  transmission  goes  to  zero,  and  finally  a  point  representative  of  conditions 
near  the  flame  front.  The  point-by-point  flame  front  position  histories  were  used  to  produce 
distances  between  the  flame  lront  and  the  fiber  at  the  time  when  the  fiber  plane  temperature  is  going 
through  these  points.  Each  of  the  points  requires  a  few  words  of  definition. 

The  definition  of  the  N20  appearance  point  is  straightforward  for  abrupt  cases  —  for 
slower  filling  cases  we  chose  the  point  at  which  the  ratio  drops  one  standard  deviation  below  its 
mean  for  the  previous  10  seconds.  A  criterion  based  on  fitted  curves  could  also  be  used,  but  those 
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Table  2.  Averages  and  Ranges  of  Quantities  Measured  in  XM39  Experiments 


1  atm  (7  cases) 

9  atm  (6  cases) 

Burning  Rate,  cm/s 

().<)  1 5()±(>.(>0 1 4 

0.041  (>±0.0025 

N-.0  Apnearance 

N20  Mole  Fraction.  Range 

(0.05-0.19) 

(0.05-0.24) 

Flame-to-Fiber  Distance  at  N^O  Appearance  Point,  cm 
From  Point-by-Point  Measurements 

From  Linear  Fit  to  Flame  Front  Position 

0.2 1  ±0. 1  1 
0.26±0. 1 1 

0.16±0. 12 
0.15±0.1 1 

Thermocouple  Readings  at  N20  Appearance.  Range,  °C 

(240-410) 

(35-155) 

Fill  Time  (from  First  N20  to  Steady  Value).  Range,  s 
(Numbers  of  Observations  in  Brackets) 

(0. 3-2.0)  (5) 

(0. 1-0.2)  [21 

(0. 1-0.3)  [51 
(0.7)  [1] 

Steady  Slate  Time  (for  Short  Fill  Time  Cases.  Time 
During  Which  a  Roughly  Constant  Transmission 
Indicating  N20  Presence  is  Maintained).  Range,  s 

(1-2) 

(0.3-0. 7) 

At,,.  Flame-to-Fiber  Distance/Buming  Rate,  s 

(8-28) 

(  1. 7-4.8) 

Change  in  a  Point  in  Temperature  Historv 

Distance  to  Flame  Front  at  a  Change  Point,  cm 

0.34±0. 17 

0.07±0.07 

Temperature  at  a  Change  Point.  °C 

I40±30 

1 10±40 

Temperature  Historv  Inflection  Point 

Distance  to  Flame  Front  at  Inflection  Point,  cm 

().26±().()85 

0.05±0.06 

Temperature  at  Temperature  History  Inflection  Point.°C 

2()()±20 

1 85±45 

Slope  at  Temperature  Inflection  Point.  °C/s 

35±1 5 

4()0±2()0 

Transmission  Loss  Point 

Flame-to-Fiber  Distance  at  Transmission  Loss,  cm 

0. 19±0.07 

().()8±().()8 

Temperature  at  Transmission  Loss  Point.  °C 

(250-500) 

(50-330) 

Expected  Time  of  Arrival  of  Flame  Front  at  Fiber  Plane 

Distance  to  Flame  Front  at  At,,  after  N20  Appears,  cm 

0.06±0.03 

-0.02±0.0(S 

Temperatures  at  At,,  after  N20  Appearance,  °C 

380±45 

410±90 

Slope  of  Temperature  History  at  At,f  after  N2Q,  °C/s 

l()±6 

100±40 
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points  (Jittered  by  less  than  1  second  front  those  derived  as  above.  As  will  be  discussed  below, 
the  rise  of  temperature  with  time  has  an  exponential  form  until,  at  an  upper  limit  temperature,  it 
departs  from  that  dependence.  We  define  that  point  as  the  temperature  history  inflection  point,  and 
defer  a  discussion  of  its  significance  to  the  following  section,  but  for  now  we  can  note  that  even 
with  the  often  uneven  burning  in  the  fiber  region,  the  average  values  of  temperature  and 
temperature  rate  of  change  at  this  point  have  fairly  small  uncertainties.  Furthermore,  in  the 
exponential  range  of  the  temperature  history,  a  discontinuous  change  in  the  parameters  of  the 
exponential  is  observed.  As  will  be  discussed  below,  this  is  most  easiliy  ascribed  to  a  physical 
change  in  the  propeliant  which  affects  the  thermal  diffusi vity  (denoted  by  the  symbol  a),  so 
Table  2  refers  to  the  a  change  point,  although  a  more  neutral  term  such  as  gradient  change  point 
might  be  preferred.  The  transmission  loss  point  is  easily  defined  as  the  point  where  one  detector 
intensity  came  so  close  to  zero  that  the  ratio  falls  well  outside  of  its  previous  range  of  variation. 

The  most  straightforward  way  to  derive  a  time  for  the  last,  flame  front  arrival  event  would 
be  to  use  the  time  at  which  the  video  flame  front  position  history  says  the  flame  front  is  at  the  fiber 
plane.  However,  we  have  already  mentioned  that  these  generally  quite  smooth  variations  of  flame 
front  positions  with  time  become  more  irregular  near  the  fiber  plane.  It  turned  out  that  evaluating 
average  quantities  at  the  zero  flamc-to-fiber  distance  time  led  to  very  wide  ranges  of  variation.  If 
we  assumed  that  the  regression  of  the  center  of  the  strand  was  more  regular,  we  could  instead  use 
the  linear  fils  to  derive  a  time  for  flame  front  arrival.  In  Table  2  we  use  a  third  method, 
extrapolating  from  a  point  where  we  expect  the  flame  front  and  burning  rate  to  be  relatively 
unperturbed.  This  point  could  be  defined  in  a  variety  of  ways,  but  we  chose  it  to  be  the  N2O 
appearance  point,  simply  because  its  flame-to-fiber  distance  was  of  intrinsic  interest  and  already 
had  been  tabulated.  This  distance  and  the  burning  rate  are  used  to  derive  a  time  between  NT) 
appearance  and  the  arrival  of  the  flame  front.  Att(,  reported  in  the  N2O  appearance  section.  This 
time  is  then  added  to  the  N2O  appearance  time  to  generate  a  time  at  which  to  evaluate  properties 
which  may  represent  conditions  near  the  flame  front.  Simply  as  a  check,  we  evaluate  the  position 
relative  to  the  flame  front  from  the  video  record  for  these  times.  As  can  be  seen  in  Table  2.  on  the 
average  the  video  flame  front  is  very  close  to  the  fiber  plane  at  At(t  from  N2O  appearance  for  the  6 
atm  cases,  but  appears  to  be  slightly  above  the  fiber  plane  in  the  1  atm  cases. 

It  is  of  interest  to  examine  the  flame-lo-liber  distance  entries  in  Table  2  to  establish  the 
sequence  of  events  experienced  at  the  fiber  insertion  plane.  The  sequence  is  different  at  the  two 
pressures.  At  1  atm.  the  temperature  inflection  point  can  arrive  first,  followed  in  fairly  rapid 
succesion  by  the  appearance  of  N2O  and  the  loss  of  transmission,  with  the  arrival  of  the  flame 
front  occurring  considerably  later.  At  6  atm,  the  N2O  appears  first,  substantially  before  the  arrival 
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o!  the  temperature  lntleetion  point,  and  therelore  in  a  lower  temperature  regime.  Loss  ol 
transmission,  inflection  point  arrival,  and  arrival  of  the  llame  ti out  region  a.s  defined  by  the  above 
procedure,  all  occur  at  about  the  same  time.  In  the  lollowing  section,  we  will  show  that  these 
changes  in  the  sequence  of  events  follow  naturally  from  the  dependence  of  the  thermal  wave  on  the 
burning  rate  and  the  assumption  that  the  length  that  NNOcan  travel  through  the  propellant  is  less 
dependent  on  the  burning  rate  and  more  a  constant  of  the  material. 


5.0  THERMAL  WAVE  ANALYSIS 


Before  we  can  hypothesize  about  the  source  of  the  NoO  we  observe,  we  must  understand 
why  it  appears  in  two  different  temperature  ranges  for  the  two  pressure  regimes.  We  can  analyze 
the  thermocouple  temperature  records  using  a  model  which  utilizes  information  on  the 
thermophysical  properties  of  the  propellant.  This  is  of  interest  not  because  it  provides  accurate 
measurements  of  these  properties  (again,  individual  bums  showed  wide  variations  in  almost  all 
observed  parameters),  but  because  it  will  be  found  that  the  thermophysical  properties  derived  from 
the  average  parameters  for  all  bums  arc  roughly  consistent  with  literature  values,  and  so  the 
expected  time/temperature  history  based  on  those  known  thermophysical  parameters  can  be  used  to 
interpret  our  results. 

It  has  long  been  known1  ^  that  solution  of  the  heat  balance  equation  for  a  thermal  wave  with 
a  constant  surface  temperature  progressing  into  a  solid  having  constant  themial  properties  and  no 
heat  sources  or  sinks  leads  to  an  expression  for  the  steady  state  temperature  profile  which  has  the 
generalized  form. 


(T-Tfl)  =  (Tinn-Tn)exp(-(x-X0)/8).  ( 1 ) 

In  this  expression,  x  is  the  (positive)  distance  into  the  propellant  from  the  flame  front.  T  is  the 
temperature  at  x.  and  Tn  is  the  ambient  temperature.  From  the  derivation.  Tmll  is  the  temperature  at 
which  the  thermal  wave  departs  from  exponential  behavior,  which  occurs  at  the  plane  where  heat 
production  (  )  becomes  important.  The  effect  of  this  heal  source  is  seen  in  the  shallower 
temperature  gradient  in  the  non-exponential  region,  a  gradient  extending  further  into  the  solid  due 
to  heat  production.  In  the  early  work  on  double-base  propellants.  Tm(1  could  be  identified  with  the 
burning  surface  temperature  Tv  It  will  be  seen  that  this  is  not  the  case  for  the  composite  propellant 
studied  here,  since  Table  2  shows  that  a  departure  from  exponential  behavior  occurs  at  a  much 
lower  temperature  than  the  surface  temperature.  This  difference  results  in  a  second  modification  to 
the  exponential  expression,  the  presence  of  the  parameter  Xft  which  is  the  distance  between  the 
flame  front  and  the  temperature  history  inflection  point.  Finally,  the  length  scale  parameter  8  is 
related  to  the  thermal  diffusivity  a  and  the  burning  rate  rh  by 

5=  aJxh  (2) 

Equation  2  represents  a  competition  between  the  ability  of  heat  to  move  through  the  solid  and  the 
speed  of  the  flame  consuming  the  solid.  If  the  former  is  large,  8  is  large  and  heat  penetrates  a  long 


IX 


distance  into  the  solid.  When  the  hunting  rate  is  relatively  large.  6  is  small  and  the  temperature 
gradient  extends  only  a  short  distance  into  the  propellant.  The  thermal  diliusivity  is  in  turn  related 
to  the  thermal  conductivity  A  and  the  heal  capacity  cp  by  a  -  A/pcp  where  p  is  the  density.  Our 
observations,  of  course,  do  not  shed  any  light  on  these  quantities  individually. 

Tire  applicability  of  a  simple  model  becomes  clearer  when  the  temperature  with  time  data 
are  displayed  in  a  semilog  plot,  as  in  Figures  10  and  1 1 .  It  is  seen  that  simple  exponential  forms  lit 
very  well  over  two  temperature  regions.  The  parameter  6  is  derived  from  such  plots  by  dividing 
the  average  bum  rate  by  the  slope  of  each  temperature  versus  lime  curve.  Therefore,  to  begin  by 
representing  the  temperature  histories  in  the  smallest  number  of  parameters,  we  fit  two  exponential 
length  scale  parameters.  8.  to  each  case,  then  computed  a  mean  and  a  standard  deviation  for  each 
pressure  and  temperature  regime.  These  parameters  are  reported  in  Table  3.  Then  we  used  the 
inflection  point  temperatures  and  positions  reported  in  Table  2  to  solve  for  a  pre -exponential  factor 
('.  m 


(T  -  T(>)  =  Cexp(-x/8)  (3) 

These  parameters  are  also  reported  in  Table  3. 

A  poin1  which  should  be  mentioned  now  is  that  the  assumption  that  the  change  in  gradient 
parameters  between  the  two  temperature  regions  corresponds  to  a  discontinuous  change  in  the 
thermal  diliusivity  (perhaps  due  simply  to  bubbles  or  cracks  appearing  at  the  transistion 
temperature)  is  only  one  possible  explanation.  Another  is  dial  the  steep  exponential  region  between 
the  two  transition  temperatures  is  due  to  an  endothermic  process  such  as  melting  or  boiling.  The 
drawback  to  this  alternative  is  that  the  exponential  form  would  not  be  expected  to  hold,  without  a 
fortuitous  dependence  of  the  rate  of  heat  loss  on  temperature. 

If  we  keep  the  hypothesis  that  the  transition  between  two  exponential  gradients  involves  a 
change  in  thermal  diliusivity.  we  must  next  address  the  fact  known  from  other  work  that  the 
thermal  diliusivity  is  a  function  of  temperature.  Experiments  in  the  -20  to  50  C  range  have  yielded 
a  three  parameter  fit  to  ouT)14  Having  shown  that  we  can  represent  our  observed  data  in  terms  of 
simple  exponential  forms,  we  now  want  to  show  that  those  data  are  consistent  with  a  plausible 
temperature  dependent  thermal  diliusivity.  For  a  temperature  dependent  a.  the  closed  form 
Equation  1  does  not  obtain.  However,  it  will  be  sufficient  to  approximate  the  true  (non-closed- 
lorm)  solution  by  Equation  I  with  a  temperature  dependent  a.  and  to  use  a  linearly  temperature 
dependent  a(T).  a(T )=u<)  +  «|T,  because  all  we  want  to  do  is  show  that  our  data  are 


Figure  10.  Temperature  History  of  6  atm  Case  Shown  in  Figure  4.  Showing  Fits  to  Simple 
Exponential  Expressions  (Dotted  Lines). 


Figure  1 1.  Temperature  History  of  1  atm  Case  Shown  in  Figure  5.  Showing  Fits  to  Simple 
Exponential  Expressions  (Dotted  Lines). 
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Table  3.  Parameters  Obtained  from  Thermocouple  Temperature  Histories 


Exponential  Fit  to  Low  Temperature  Range 

1  utm 

6  atm 

8.  Length  Scale  Parameter,  cm 

0.165±0.()32 

().()42±0.0065 

C.  Pre-exponential  Factor  (Uncertainty  Range ).°C 

595  (430-620) 

330  (185-430) 

Exponential  Fit  to  Hieh  Temperature  Range 

8.  Length  Scale  Parameter,  cm 

0.095±0.044 

0.()24±0.009 

C,  Pre-exponential  Factor  (Uncertainly  Range).  °C 

2680(1260-25200) 

1090  (825-2470) 

Model  Parameters  Consistent  with  Exponential  Fit 

Ti„,rT().  T 

174 

160 

Thermal  Wave  Origin.  Xn,  cm 

0.26 

0.045 

Low  Temperature  Range 

Thermal  DilTusivity  Constant  Term.  a<),  cm*Vs 

0.0019 

a  Temperature  CoelTicient.  exj .  cm2/s  "O' 

9  x  KK1 

High  Temperature  Range 

Thermal  DilTusivity  Constant  Term.  a<>  cm2/s 

0.0014 

a  Temperature  CoelTicient.  a|.  em’/s  °C'1 

2  x  l()-fl 
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noisy  enough  that  it  is  difficult  to  distinguish  between  a  constant  and  a  temperature-dependent 
thermal  dil'fusivity. 

To  do  this,  we  use  the  Equation  3  representations  of  average  temperature  histories,  together 
with  upper  and  lower  error  limit  expressions  derived  as  follows.  The  6  parameter  in  the  upper 
limit  expression  is  the  average  8  plus  its  standard  deviation.  The  C  parameter  for  the  high 
temperature  range  is  then  obtained  by  fitting  to  the  inflection  point  temperature  plus  its  standard 
deviation.  The  C  parameter  for  the  low  temperature  range  is  set  by  matching  the  high  temperature 
expression  at  the  a  change  point.  Since  the  8  parameters,  which  appear  in  the  exponential,  have 
fairly  large  standard  deviations,  and  since  the  range  of  observed  temperatures  at  any  distance  from 
the  flame  front  is  relatively  small,  the  C  parameters  are  forced  into  some  large  variations,  as  seen  in 
the  ranges  reported  in  Table  3. 

To  show  how  these  exponential  forms  compare  with  some  individual  observations. 

Figures  12  and  13  once  again  exhibit  the  temperature  data  for  the  two  cases  shown  in  Figures  10 
and  1 1 .  Here,  however,  the  video  flame  front  position  record  has  been  used  to  convert  the  time 
axis  into  a  distance  from  flame  front  axis.  The  exponential  data  representations  using  the 
parameters  in  Table  3,  and  the  uncertainty  limit  expressions  derived  as  discussed  above,  are  also 
plotted  in  the  figures.  It  can  be  seen  that  the  substantial  variations  in  the  distance  scales  cause  these 
observed  curves  to  fall  outside  the  error  limit  expressions.  In  the  6  atm  case,  the  observed 
inflection  and  a  change  points  occur  at  0. 10  and  0. 1 3  cm  from  the  flame  front,  much  further  out 
than  the  average  values  of  0.05  and  0.07  cm.  In  the  1  atm  case,  on  the  other  hand,  inflection  and 
a  change  point  at  0.22  and  0.26  cm  are  well  inside  the  average  values  of  0.26  and  0.34  cm. 
However,  the  uncertainty  in  the  position  axis  does  not  directly  affect  the  determination  of  thermal 
dil'fusivity  values  from  the  slope  of  the  temperature  history  curve. 

Now  that  we  have  a  compact  way  of  expressing  what,  on  average,  the  temperature  histories 
were  in  the  two  pressure  regimes,  we  can  go  on  to  compare  to  expected  temperature  histories  as 
predicted  by  Equation  1.  In  Figures  14  and  15  we  present  the  Equation  3  exponential  data 
representations  as  dashed  lines,  with  the  error  limits  expressions  again  yielding  two  bracketing 
dotted  lines.  The  solid  lines  are  given  by  Equation  1  using  the  parameters  in  the  lower  half  of 
Table  3.  The  pre-exponential  factor  C  can  be  divided  into  the  inflection  point  temperatures  given  in 
Table  2  and  exponentials  involving  Xo  parameters  which  are  essentially  the  average  observed 
positions  of  those  inflection  points.  The  linear  a(T)  expressions  given  by  the  parameters  in 
Table  3  are  those  which  minimize  the  squared  differences  between  the  dashed  and  solid  lines  in 
Figures  14  and  15.  These  linear  thermal  diffusivity  expressions  are  plotted  in  Figure  16.  along 
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Figure  12  Thermocouple  Reading  versus  Flame  Front  Position  of  6  atm  Case  Shown  in  Figure  4 
(Solid  Line),  Average  Exponential  Data  Representations  (Dashed  Lines),  and  Error 
Limit  Expressions  (Dotted  Lines). 


Figure  13.  Thermocouple  Reading  versus  Flame  Front  Position  of  1  atm  Case  Shown  in  Figure  5. 
(Line  types  have  same  definitions  as  in  Figure  13.) 
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Figure  14.  Comparing  Simple  Exponential  Data  Representation  and  Error  Limit  Expressions  for 
6  atm  Cases  (Dashed  and  Dotted  Lines)  with  Equation  1  Expression  Including 
Temperature-Dependent  Thermal  Diffusivity  and  Other  Table  3  Parameters  (Solid 
Lines). 


Figure  15.  Comparing  Simple  Exponential  Data  Representation  and  Error  Limit  Expressions  for 
1  atm  Cases  (Dashed  and  Dotted  Lines)  with  Equation  1  Expression  Including 
Temperature-Dependent  Thermal  Diffusivity  and  Other  Table  3  Parameters  (Solid 
Lines). 
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Figure  16.  Thermal  Dil tusivity  tor  XM39:  Low  Temperature  Measurements14  (Solid  Line).  Data 
Point  from  Strand  Burning  Observations  of  Reference  15  (With  Quoted  Error  Limits, 
Plotted  at  Observed  Inflection  Point  Temperature  of  165  °C),  and  Our  Linear  a(T) 
Expressions  From  Table  3,  Used  in  Figures  14  and  15.  (Dot-Dash  and  Dashed  Lines). 

parameters  which  give  the  best  match  between  Equation  1  and  Equation  3  expressions  while  still 
being  roughly  compatible  with  the  other  measurements  of  a.  many  other  a(T)  values  would  also 
fit  our  observations  well  enough.  This  is  shown  in  Figure  17.  which  shows  a  region  of  a(T) 
values  which  give  temperature  history  curves  which  lie  within  the  error  limit  expressions  for  both 
the  1  atm  and  6  atm  cases.  In  other  words,  the  uncertainty  in  what  values  of  a(T)  are  consistent 
with  our  observations  is  at  least  as  large  as  indicated  in  Figure  17.  Considering  these  large 
uncertainties,  it  is  not  a  matter  of  great  concern  that  our  average  a  values  at  low  temperatures  seem 
to  be  higher  than  those  of  Reference  9.  nor  that  our  high  temperature  a  values  seem  to  be  below 
that  reported  in  Reference  15. 

To  summarize,  a  single  set  of  a(T)  expressions  can  be  used  to  reproduce  the  average 
temperature  histories  for  both  the  1  atm  and  6  atm  data  sets.  This  is  in  spite  of  the  fact  that  (at  the 
simple  exponential  level  of  data  analysis)  the  average  6  parameters  for  the  two  data  sets  are  not  in 
the  ratio  of  the  burning  rales  for  the  two  pressures  (rh(6  atm)/rh(  1  atm)  =  2.7  ±  0.3),  but  instead, 
the  ratios  in  both  temperature  regimes  are  about  4.  However,  the  uncertainties  in  the  ratios  of 
6(6  atm)  to  5(  1  atm)  are  large  (±1.0  and  2.4  for  low  and  high  temperature  ranges),  so  much  of  the 
difference  can  be  put  down  to  large  standard  deviations,  with  a  component  also  potentially 
ascribablc  to  the  temperature  dependence  of  a  and  the  fact  that  somewhat  different  ranges  of 
temperature  contribute  to  single  values  of  a  derived  from  the  two  different  pressure  cases. 
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To  summarize,  a  single  set  of  a(T)  expressions  can  he  used  to  reproduce  the  average 
temperature  histories  for  both  the  1  atm  and  6  atm  data  sets.  This  is  in  spite  of  the  fact  that  (at  the 
simple  exponential  level  of  data  analysis)  the  average  8  parameters  for  the  two  data  sets  are  not  in 
the  ratio  of  the  burning  rates  for  the  two  pressures  (rh(6  atm)/rb(  1  atm)  =  2.7  ±  0.3),  but  instead, 
the  ratios  in  both  temperature  regimes  are  about  4.  However,  the  uncertainties  in  the  ratios  of 
8(6  atm)  to  8(  1  atm)  arc  large  (±1.0  and  2.4  for  low  and  high  temperature  ranges),  so  much  of  the 
difference  can  be  put  down  to  large  standard  deviations,  with  a  component  also  potentially 
ascribablc  to  the  temperature  dependence  of  a  and  the  fact  that  somewhat  different  ranges  of 
temperature  contribute  to  single  values  of  a  derived  from  the  two  different  pressure  cases. 
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Figure  17.  Thermal  Diflusivity  for  XM39:  Low  Temperature  Measurements14  (Solid  Line).  Data 
Point  from  Strand  Burning  Observations1- ,  and  Range  of  Values  whose  Equation  1 
Temperature  Histories  Fall  within  the  Error  Limit  Regions  in  Figures  14  and  15  (Dots). 
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6.0  DISCUSSION 


This  study  has  produced  three  major  results.  The  first  is  that  we  have  observed  N20.  a 
easeous  decomposition  product  of  RDX.  evolving  into  a  small  volume  in  the  condensed  phase  of  a 
burning  nitramine  propellant  strand.  As  we  will  discuss  below,  this  result  is  expected,  in  that  it  is 
consistent  with  other  observations  of  the  chemistry  of  RDX  decomposition.  The  second  is  that  in 
the  1  atm  cases,  there  is  strong  evidence  that  N20  is  not  seen  until  substantially  alter  the  RDX 
melting  and  decomposition  processes  have  begun.  It  has  been  expected  that  N20  is  not  a  primary 
product  of  RDX  decomposition,  nut  it  is  somewhat  surprising  to  see  such  a  clear  separation.  The 
third  result  is  that  many  6  atm  cases  show  an  abrupt  N20  evolution  at  a  time  when  the  propellant 
forming  the  observation  volume  is  still  at  a  temperature  well  below  that  needed  to  produce  rapid 
decomposition  of  RDX.  We  can  suggest  a  mechanism  for  this  observation,  but  first  we  need  to 
place  it  in  the  context  of  what  we  and  others  have  learned  about  XM39  burning  and  RDX 
decomposition. 

6. 1  Correlation  of  Observed  Events  with  Temperature  History 

To  aid  in  visualizing  the  relationships  between  the  events  listed  in  Table  2.  in  Figure  18  we 
present  a  sketch  of  temperature  histories  and  event  positions  and  temperatures  which  arc  roughly 
consistent  with  those  in  the  table.  It  turns  out  that  it  is  not  possible  to  construct  a  self-consistent 
graph  out  of  the  averages  listed  in  Table  2.  because  the  large  variances  and  small  numbers  of  points 
mean  that  different  runs  have  different  weights  in  determining  the  various  parameters.  Even  such 
basic  qualitative  characteristics  as  the  order  of  events  cannot  be  completely  summarized  by  a  single 
figure.  For  example,  although  four  6  atm  runs  do  indeed  have  the  order  shown  in  the  figure. 

N20  appearance/transniission  loss/a  change/inflection  point,  two  others  have  a  different  order, 
a  change/N20  appcarance/inflection  point/transmission  loss.  Still,  to  the  extent  that  there  is  a 
typical  set  of  positions  and  temperatures  associated  with  the  set  of  observed  events.  Figure  18  is  as 
good  a  representation  of  it  as  our  present  data  set  will  allow. 

In  the  context  of  Figure  18,  then,  we  want  to  ask  what  we  know  about  the  state  of  the 
propellant  through  the  time  interval  surrounding  the  observation  of  N20  appearance  in  the 
observation  region,  based  both  on  our  data  and  on  other  sources.  We  begin  by  noting  that  in 
Table  2.  three  of  the  events  are  associated  with  temperatures  with  fairly  small  standard  deviations, 
suggesting  that  they  involve  processes  which  are  legitimately  associated  with  specific  temperatures. 
The  burning  surface,  of  course,  is  expected  to  have  a  temperature  which  is  only  secondarily  a 
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DISTANCE  FROM  FLAME  FRONT,  cm 

Figure  1 X.  Qualitative  Sketch  of  Typical  Thermal  Waves  and  Positions  of  Other  Events  in 
Strand  Burning  Observations 

function  of  combustion  pressure  (although  our  strand  burning  in  air  is  certainly  expected  to  have  a 
higher  surface  temperature  than  observed  in  inert  atmospheres).  The  value  from  our  extrapolation 
procedure  for  the  time  of  arrival  of  the  flame  front,  of  about  400  °C.  is  consistent  with  the  range 
quoted  by  Brill  and  Brush16  of  350-400  °C  for  the  surface  temperature  of  burning  HMX  and  RDX. 

Again,  the  temperature  inflection  point,  occurring  around  190  °C.  is  thus  not  to  be 
identified  with  the  burning  surface  temperature.  This  conclusion  is  reinforced  by  the  fact  that  in  1 
atm  cases  we  have  good  infrared  transmission  over  a  long  path  through  a  channel  in  the  solid  long 
after  the  inflection  point  has  passed  the  observation  region.  We  know  the  major  component.  RDX. 
does  not  decompose  at  significant  rates  until  it  has  melted.  Consulting  Table  1.  we  find  that  RDX 
melts  at  190  °C.  (This  value  is  a  manufacturer's  specification — the  references  to  values  in  the  203- 
2()5°C  range  may  refer  to  material  of  higher  purity  than  would  be  found  in  gun  propellant-2-617.) 
Therefore,  our  inflection  point  may  be  associated  with  the  melting,  decomposition,  and  eventual 
heat  releasing  processes  in  RDX.  It  is  clear  that  if  RDX  melting  begins  at  the  inflection  point 
temperature,  its  effect  on  the  temperature  history  is  masked  by  the  exothermic  process.  An 
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alternative  explanation,  that  what  we  have  referred  to  as  the  a  change  point  is  actually  the  RDX 
melting  point,  would  shift  all  our  observed  temperatures  to  surprisingly  high  values. 


Because  these  processes  take  place  at  finite  rates,  as  the  burning  rate  increases  we  expect 
the  inflection  point  to  shift  to  higher  temperatures,  as  seen  in  the  observations  of  Reference  15. 
These  covered  the  pressure  range  from  75  to  750  psig.  and  showed  a  strong  dependence  of  the 
inflection  point  on  pressure  The  temperatures  reported  at  the  two  extreme  pressures  are  about  165 
and  370  °C.  respectively,  while  at  2(X)  and  5(X)  psig  the  temperatures  are  250  and  300  °C.  It  can  be 
seen  that  the  expected  variation  from  1  to  6  atm  is  smaller  than  our  standard  deviations,  so  our 
failure  to  observe  this  trend  is  not  surprising. 

Finally,  we  want  to  ask  if  the  lower  temperature  change  in  the  form  of  the  thermal  wave, 
which  we  have  referred  to  as  the  a  change  point,  can  be  identified  with  a  property  of  one  of  the 
propellant  ingredients.  If  we  note  the  temperature  ranges  in  Table  2.  centered  around  1 25  °C.  and 
then  consult  Table  1 ,  we  find  the  closest  temperature  listed  to  be  the  boiling  point  of  the  plasticizer. 
Acetyl  triethyl  citrate  is  a  liquid  at  room  temperature,  and  its  plasticizing  action  essentially  involves 
partial  solvation  of  the  other  propellant  ingredients.  The  thermal  properties  of  this  solution  will  be 
different  than  those  of  the  pure  liquid  plasticizer,  but  it  is  not  unreasonable  to  expect  a  sharp  change 
in  properties  near  the  plasticizer  boiling  point,  due  to  either  vaporization  or  thermal  decomposition. 
Here  again  it  must  be  remarked  that  in  their  observations  of  strand  burning  of  XM39  with 
embedded  thermocouples,  neither  Reference  9  or  15  noted  such  a  point  in  the  temperature  history. 
We  can  only  suggest  that  at  their  higher  pressures  and  faster  time  scales,  the  process  involved  is 
too  slow  to  be  observed  as  a  sharp  discontinuity.  This  point  too  is  not  seen  by  us  to  move  towards 
higher  temperatures  at  faster  burning  rates,  but  again  our  error  limits  are  substantial.  We  should 
also  note  again  that  the  steeper  gradient  between  inflection  and  a  change  points  may  not  be  a 
thermal  diflusivity  effect  at  all.  but  simply  reflect  a  heat  consumption  process  such  as  melting, 
boiling,  or  thermal  decomposition. 

Having  made  the  case  that  three  of  the  events  identified  in  Table  2  and  Figure  IS  are,  under 
our  conditions,  correlated  with  relatively  fixed  temperatures,  what  can  be  said  of  the  other  two 
events?  It  is  clear  that  transmission  loss  is  not  correlated  with  temperature,  and  the  wide  variation 
in  temperatures  at  which  it  occurs  tells  us  that  a  number  of  different  processes  can  be  involved  in 
bringing  the  transmission  to  zero,  including  misalignment  of  the  fibers  and  intrusion  of  molten 
propellant  material  into  the  observation  region.  It  turns  out.  however,  to  be  fairly  well  correlated 
with  distance  from  the  flame  front,  with  the  standard  deviation  in  position  for  the  1  atm  cases  being 
particularly  small. 
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Finally,  u  is  clear  from  Figure  IS  that  appearance  ol  NAO  in  the  observation  volume  is  not 
correlated  with  the  temperature  at  that  point.  If  anything,  it  too  is  better  correlated  with  distance 
trom  the  flame  front.  To  understand  why  we  suggest  that  the  NAO  we  observe  has  been  formed 
not  in  the  observation  region,  but  substantially  nearer  the  tlame  front,  we  need  to  review  what  we 
know  about  combustion  of  nitrannne  propellants  like  XN139. 

First,  we  know  that  NiO  is  a  product  of  thermal  decomposition  of  RDX.  but  not  a  final 
product  of  the  gas  phase  combustion  of  those  solid  decomposition  produces,  either  in  the  primary 
tlame  or  in  secondary  combustion  with  added  oxygen.  The  evolution  of  N^O  from  heated  RDX 
has  been  observed  by  a  number  of  investigators,  using  mass  spectrometric  and  infrared  absorption 
detection  techniques.-^61 116-1X  An  example  of  the  quantitative  measurements  of  Fetherolf  and 
Litzinger 1 1  is  given  in  Figure  19.  along  with  our  estimate  of  how  the  somewhat  broader  range  of 
species  observed  by  Behrens1*  in  his  low  pressure  experiments  would  evolve  under  the  I  atm 
conditions  of  Reference  1 1 .  The  point  of  this  figure  is  that  under  a  wide  range  of  conditions  N->0 
is  expected  to  be  a  major  gas  phase  product  of  the  thermal  decomposition  of  RDX.  The  second 
point,  that  NoO  is  not  expected  to  be  a  major  product  of  gas  phase  combustion  of  RDX 
decomposition  products,  is  made  by  Figure  20.  It  shows  results  of  the  gas  phase  chemistry 
calculations  of  Melius|y-20  for  RDX  burning  in  I  atmosphere  of  an  inert  gas.  In  his  mechanism. 
NSO  is  consumed  by  thermal  decomposition  as  well  as  reaction  with  hydrogen  atoms,  and  has 
essentially  disappeared  0. 1  cm  above  the  burning  surface.  Thus,  we  know  that  the  gas  we  are 
observing  comes  from  inside  the  condensed  phase  of  the  propellant,  not  from  the  exhaust  gas  from 
the  flame.  In  any  case,  the  rapid  gas  tlow  over  the  strand  prevents  recirculation  of  the  exhaust  gas 
into  the  1 1  her  insertion  region. 

We  also  know  that  measurements  of  the  kinetics  of  RDX  decomposition  show  that  gas 
evolution  is  very  slow  below  the  melting  point. •s-lx--1  Yet  we  also  know  that  infrared  transmission 
through  even  0. 1  cm  of  molten  propellant  will  be  essentially  zero,  so  that  the  walls  of  the 
observa'ion  region  must  be,  at  least,  below  the  melting  point  of  the  CAB  binder  (215-225  T). 
Further,  we  know  that  the  NiO  appearance  rale  in  cases  like  those  of  Figures  4  and  5  is  faster,  by 
one  to  two  orders  of  magnitude,  than  could  be  expected  from  NAO  evolution  from  the  walls  of  the 
observation  region,  when  the  observed  evolution  rate  from  Reference  IS  for  RDX  at  205  C  is 
used.  All  these  facts  lead  us  to  suggest  that  NSO  from  RDX  decomposition  accumulates  in 
pressurized  bubbles  in  regions  of  the  propellant  strand  which  are  hotter  and  closer  to  the  flame 
tront  than  the  observation  region,  and  that  the  NiO  is  transported  from  the  pressurized  bubbles  to 
the  observation  volume  by  the  opening  of  cracks  connecting  the  two  regions. 


M) 


Figure  19. 


Major  RDX  Decomposition  Gas  Species  from  Experiments  of  Fetherolf  and 
Lit/.inger. 1 1  Compared  to  Adaptation  of  Observations  of  Behrens  (dots)lx 


Figure  20. 


Major  Gas  Phase  Products  of  RDX  Primary  Combustion,  from  Model  Calculations 
of  Melius.'1^'0 


Several  groups  have  made  scanning  electron  microscope  <SEM)  observations  ot  sections 
through  the  surface  oi  quenched  strands  ot  XM.W --  A  melt  region  is  observed,  ot  thickness 
(MM -0.03  cm  lor  burning  in  air  at  1  atmosphere  pressure.  This  melt  region  is  tilled  with  bubbles 
ot  various  si/es.  ranging  in  diameter  up  to  the  lull  thickness  ot  the  melt  layer.  Because  these 
bubbles  are  so  close  to  the  burning  surface,  it  would  be  surprising  it  they  were  the  source  of  NAO 
observed  on  the  order  of  0.1  cm  further  down,  in  relatively  still-solid  propellant.  We  would  expect 
most  of  these  bubbles  to  burst  upward,  ejecting  decomposition  gas  into  the  flame,  and  indeed  what 
we  could  see  of  our  burning  surfaces  seemed  to  be  consistent  with  a  melt  layer  containing  growing 
and  bursting  bubbles.  Furthermore,  a  single  bubble,  even  with  a  0.03  cm  (300  pm )  diameter,  will 
not  till  the  observation  region,  a  cylinder  of  0.02  cm  diameter  and  order  of  0.3  cm  length,  to  the 
mole  fraction  levels  we  observe,  which  are  consistent  with  pure  decomposition  gas  at  ambient  ( 1  or 
b  atm)  pressure. 

These  objections  lead  to  the  suggestions  that  the  bubbles,  or  pressurized  voids,  which  are 
the  source  of  the  NAO  we  observe,  are  not  in  the  melt  region  but  below  it.  in  a  region  where  at  least 
the  higher-melting  binder  has  kept  enough  strength  to  hold  these  bubbles  in  place,  and  that  our 
observations  of  sudden  filling  of  die  observation  volume  with  NAO  involves  the  nearly 
simultaneous  opening  of  several  cracks  connecting  it  to  several  bubbles.  Such  cracks  have  not.  to 
our  knowledge,  been  observed  in  any  SEM  photographs  of  XM31)  strands,  but  there  need  not  be 
many  of  them,  and  they  will  be  more  difficult  to  see  than  the  bubbles.  Finally,  we  cannot  at  this 
point  rule  out  the  possibility  than  our  fiber  insertion  procedure,  involving  pressing  a  groove  at  least 
•  MM  cm  deep  into  the  propellant,  has  generated  fractures  extending  many  times  that  depth  into  the 
solid,  or  has  at  least  set  up  strains  which  make  crack  formation  more  likely  during  burning. 

o  2  Correlation  of  Our  NAO  Observations  with  Chemical  Mechanisms  ot  RDX  Decomposition 

This  lost  part  of  our  discussion  has  to  do  with  the  status  of  NAO  as  a  major  decomposition 
product  ot  RDX.  It  turns  out  that  the  observation  of  substantial  NAO  is  completely  consistent  with 
the  existing  literature  of  experiments,  theoretical  predictions,  and  analysis.  To  show  this,  we 
present  a  briel  review  of  some  of  the  major  contributions  to  the  current  overall  picture  of  RDX 
decomposition. 

It  has  been  w  idely  noted1  lx-  that  there  appear  to  be  two  global  decomposition  channels 
lor  RDX.  one  leading  to  the  products  NAO  +  IDCO.  and  one  leading  to  N(F  +  IDCN  (or 
MONO  +  HCN.  or  NAM  +  HCN  +  H.  the  differences  expressing  the  tact  that  some  of  the  products 
ot  this  channel  are  less  stable).  The  pathways  to  these  products  can  include  a  number  of  parallel 
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processes,  some  of  which  can  involve  a  number  of  sequential  steps.  Not  surprisingly,  then,  some 
variation  in  relative  amounts  ol'  the  above  products  and  others  has  been  observed  depending  on  the 
conditions  during  the  decomposition.  Experiments  on  RDX  decomposition  under  relatively 
controlled  conditions,  including  the  use  ot  deuterium  isotope  techniques,  have  unraveled  some  of 
the  details  of  these  processes. 

Behrens  and  Bulusu 1X  have  been  able  to  explain  their  simultaneous  thermogravimetry 
modulated  beam  mass  spectrometry  observations  of  RDX  thermal  decomposition  under  high 
vacuum  conditions  using  four  pathways: 

1 .  Formation  of  l-nitroso-3.5-dinitrohexahydro-.v-triazine  (ONDNTA  —  RDX  with 
one  NOt  group  replaced  by  NO)  followed  by  ring  breaking  to  yield  N^O.  CH2O, 
and  other  produces. 

2.  MONO  elimination  to  form  oxy-.v-tria/.ine  (OST).  a  molecule  which  retains  the  C-N 
ring  of  RDX  and  has  the  formula  C^H^N^O  and  an  uncertain  structure,  as  well  as 
H20.  NO.  and  NO:. 

3.  N-NO>  bond  breaking  leading  to  ring  breaking  and  formation  of  NO->  +  H-.CN  + 
2N2O  +  2CH2O. 

4.  Catalytic  decomposition  to  form  N2O.  CH20.  NCh  and  NFhCHO. 

They  find  no  evidence  that  process  4  is  the  water  catalysis  suggested  by  Melius-^,  and  instead 
suggest  the  catalyst  is  a  residue  of  RDX  decomposition  which  builds  up  on  the  walls  of  their 
reaction  cell.  That  makes  this  mechanism  specific  to  their  apparatus,  although  analogous  processes 
might  occur  in  the  burning  of  an  RDX  composite  propellant.  The  other  three  mechanisms  all 
involve  N-N  bond  breaking,  since  process  1  is  thought  to  proceed  by  N-NCA  bond  cleavage 
followed  by  reaction  with  a  source  of  NO.-y  This  might  seem  to  be  in  conflict  with  the  deuterium 
kinetic  isotope  el  lect  studies  of  Shackelford  et «/. 30.3 1  who  found  C-H  bond  rupture  to  be  the  rate- 
controlling  step.  However.  Melius-X  has  pointed  out  that  not  only  does  N-NCH  bond  breaking 
weaken  remaining  bonds,  but  C-H  bond  breaking  results  in  a  negative  N-NCH  bond  energy  with 
only  a  2  kcal-mole*1  barrier  for  leaving.  Furthermore.  Behrens  and  Bulusu  found  evidence  of  C-H 
bond  breaking  in  both  pathways  1  and  2.  which  together  result  in  two-thirds  of  the  total  RDX 
decomposition  in  their  experiment.  Behrens  and  Bulusu  also  note  that  pathway  3.  whose 
contribution  is  estimated  to  be  only  1 (>'/?•  of  the  total,  may  be  the  concerted  symmetric  triple  fission 
to  form  three  methylenenitramine  fragments  observed  in  the  gas  phase  by  Zhao  ct  <,l22.  or  may 
instead  be  a  process  yielding  two  methylenenitramine  molecules  plus  H2CN  and  NCK 


To  summarize,  the  above  experimental  evidence  seems  to  give  a  coherent  picture  which 
agrees  wuh  the  prediction  of  Melius2X  on  energetic  grounds  that  RDX  decomposition  in  the 
condensed  phase  will  be  dominated  by  processes  involving  N-NOi  bond  breaking,  at  least  at  high 
heating  rates  where  catalytic  processes  are  not  important.  The  dependence  of  processes  and 
product  distributions  on  heating  rates  is  part  of  a  more  general  dependence  on  the  environment  of 
the  RDX  molecule,  including  not  only  the  temperatures  and  time  scales  involved  but  also  the  phase 
(process  1  is  thought  to  be  dominant  in  the  solid  phase,  and  at  low  temperatures  in  general). 

Table  4  gives  a  review  of  experimental  observations  of  RDX  decomposition  to  give  NAO  ordered 
hv  heating  rate.  It  can  be  seen  that  the  only  experiments  in  which  conditions  did  not  allow 
progress  beyond  the  initial  NCb  loss  step  are  those  of  Wight  and  Botcher2-'5'27  (although  with 
longer  laser  pulses  they  could  observe  formation  of  other  products  including  NoO).  As  to  the  rest 
of  the  experiments,  covering  five  orders  of  magnitude  in  heating  rate,  it  is  not  easy  to  discern  any 
systematic  trends  in  product  distributions,  given  that  differences  in  observed  concentrations  are 
only  at  the  factor  of  two  level  and  that  other  differences  in  the  experimental  conditions  can  also 
easily  generate  differences  at  the  same  level.  Our  experiments  (with  heating  rates  taken  from  the 
temperature  history  slopes  in  Table  2)  are  seen  to  be  in  the  middle  of  this  range  of  heating  rates. 

We  can  add  the  comment  that  the  nature  of  the  initiating  process  for  RDX  thermal 
decomposition  has  been  particularly  controversial  because  of  a  seeming  disagreement  between  the 
above-mentioned  results  of  infrared  multiphoton  dissociation  experiments  by  Zhao  a  air  2 
showing  that  two-thirds  of  their  RDX  decomposed  by  a  concerted  detrimerization  process,  and  the 
condensed  phase  experiments  which  find  little  or  no  evidence  for  such  a  process.  However,  the 
suggestion  has  been  made  by  Wight  and  Botcher27  that  in  tact  no  discrepancy  exists,  that  while  the 
concerted  process  proceeds  easily  in  the  vapor  phase,  the  condensed  phase  constrains  the  motion 
of  the  RDX  molecule  too  much  for  it  to  contribute  significantly.  Thus  the  picture  of  RDX 
decomposition  in  the  condensed  phase  remains  one  in  which  all  ring-breaking  proceeds  by  several 
sequential  steps  and  results  from  the  destabilizing  effects  of  NCH  loss. 

Our  experimental  observations  of  NoO  are  consistent  w'ith  the  above  picture.  Our  heating 
rates  are  slow  compared  to  chemical  reaction  rates,  and  our  NiO  originates  in  a  mixed  gas/liquid 
region  w ith  plenty  of  opportunity  for  homogeneous  or  heterogeneous  reactions  to  carry  initial 
products  of  decomposition  to  secondary  products.  A  decomposition  product  whose  observation 
would  tell  us  more  about  the  environment  in  the  decomposition  region  and  during  their  transport 
through  the  propellant  is  NCK  which  is  expected  to  be  more  reactive  as  well  as  being  formed  by 
some  of  the  first  bond-breaking  processes  in  RDX.  Therefore,  in  ongoing  work,  we  plan  to  look 
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Brill  and  Brush  T-Jump  FTIR 

Behrens  and  Bulusu  Thermogravimetry/MBMS 
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parameters  remaining  much  the 
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Ref. 

Heating  Rate.  K/s 

NoDObsc 

2 

2.5 

Yes 

3 

8-200 

Yes 

23 

100-350 

Yes 

16 

<  2000 

Yes 

18 

0.017 

Yes 

24 

0.17 

Yes 

This 

10-100 

Yes 

Work 

25 

1-3  x  107 

No 

7.0  SUMMARY 


We  draw  three  major  conclusions  from  this  work:  NoO  is  observed  as  a  decomposition 
product  of  RDX  under  actual  propellant  combustion  conditions,  it  is  formed  in  high  temperature 
regions  of  the  condensed  phase  and  is  able  to  penetrate  into  regions  of  the  propellant  which  are 
otherwise  little  affected  by  the  advancing  llame  front,  and  under  our  conditions  its  formation  is  not 
an  initial  step  in  RDX  decomposition.  Of  course  our  observations  by  themselves  will  not  fully 
specify  the  structure  of  a  burning  strand  and  the  physical  mechanisms  which  lie  behind  it,  but  they 
certainly  strongly  encourage  the  speculation  that  the  suddenly  appearing  N2O  originates  in 
pressurized  bubbles  inside  liquid  RDX  which  in  turn  is  contained  inside  a  matrix  formed  by  the 
propellant  binder.  This  matrix  would  have  to  be  structurally  sound  enough  that  some  bubbles  are 
prevented  from  releasing  their  contents  upward  into  the  llame  region,  but  instead  eject 
decomposition  gas  into  the  observation  volume  when  cracks  open  creating  vent  paths. 
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